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ABSTRACT 

We observed NGC 1624-2, the O-type star with the largest known magnetic field 
(B p ~ 20 kG), in X-rays with the ACIS-S camera onboard the Chandra X-ray Ob¬ 
servatory. Our two observations were obtained at the minimum and maximum of the 
periodic Hct emission cycle, corresponding to the rotational phases where the magnetic 
field is the closest to equator-on and pole-on, respectively. With these observations, 
we aim to characterise the star’s magnetosphere via the X-ray emission produced by 
magnetically confined wind shocks. Our main findings are: (i) The observed spectrum 
of NGC 1624-2 is hard, similar to the magnetic O-type star 9 1 OriC, with only a few 
photons detected below 0.8 keV. The emergent X-ray flux is 30% lower at the Ha 
minimum phase, (ii) Our modelling indicated that this seemingly hard spectrum is in 
fact a consequence of relatively soft intrinsic emission, similar to other magnetic Of?p 
stars, combined with a large amount of local absorption (~l-3x 10 22 cm -2 ). This com¬ 
bination is necessary to reproduce both the prominent Mg and Si spectral features, 
and the lack of flux at low energies. NGC 1624-2 is intrinsically luminous in X-rays 
(log L ~33.4) but 70-95% of the X-ray emission produced by magnetically confined 
wind shocks is absorbed before it escapes the magnetosphere (logLx Mc ° r ~32.5). (iii) 
The high X-ray luminosity, its variation with stellar rotation, and its large attenuation 
are all consistent with a large dynamical magnetosphere with magnetically confined 
wind shocks. 

Key words: stars: magnetic fields - stars: individual: NGC 1624-2 - stars: massive 
- stars: winds, outflows - X-rays: stars - stars: early-type. 


1 INTRODUCTION 


A small fraction of O stars have large-scale, stable surface 
magnetic fields, as revealed in only the past decade by spec- 


tropolarimetric measurements of the Zeeman effect (e.g. Do- 


nati et al. 2006 Martins et al. 2010 Wade et al. 2012bI. 


Even before these direct detections, an oblique dipole field 
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was suspected in the prototype magnetic O star, 0 1 OriC, 
notably due to its luminous, rotationally modulated X-ray 
emission (Stahl et al. 1996 Gagne et al. 19971. Indeed, 
strong and rotationally modulated X-ray emission due to 
a combination of occultation and local absorption was pre¬ 
dicted to be the hallmark of the wind-fed dynamical magne¬ 
tospheres of magnetic O stars with dipole fields misaligned 


with their rotational axis (e.g. Babel & Montmerle 1997b 


Gagne et al. 2005 ud-Doula et al. 2014). The X-rays are 


due to the radiative cooling of plasma that is shock-heated 
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Figure 1. Schematic of a magnetic massive star dynamical mag¬ 
netosphere (e.g. |Sundqvist et ah]|2012| |Petit et al.||2013) . Solid 
blue lines indicate regions below the last closed magnetic loop 
that confine the wind, located near the Alfven radius R ,\. Most 
of the Ha emission originates here. Dashed lines indicate regions 
where the momentum of the wind results in open field lines. The 
bulk of the X-rays are produced in the region indicated in purple; 
see f|6] The insets illustrate the view of an observer as the star’s 
rotation changes the orientation of the magnetosphere. It is im¬ 
portant to note that due to the long rotation periods of magnetic 
O-type stars, the dynamical effects of rotation on the magneto- 
spheric structure are negligible (|ud-DouIa et al.|2008|l. 


when wind-flows from opposite hemispheres collide near the 
magnetic equatorial plane (see Fig.[T|) . Measurements of this 
X-ray emission are therefore a primary diagnostic of the 
physical conditions in, and processes that control, the mag¬ 
netospheres of 0 stars. 

The rotational modulation of the X-ray emission in 
8 1 Ori C is due to occultation of the post shock material by 
the stellar disk during equator-on viewing phases, with no 
conclusive evidence in phase-resolved high-resolution Chan¬ 
dra spectra for any large absorption of X-rays by the mag¬ 
netosphere, despite the fact that absorption by cooled post¬ 
shock plasma was expected (|Babel fe Montmerle||1997a| 
Gagne et al. 2005 |Donati et al. 20021. The surface field 


strength of a little over 1 kG implies a maximum closure ra¬ 
dius of magnetic field linetQof only slightly more than 2R*. 
The modest size of this star’s magnetosphere therefore leads 
to relatively low absorption column densities. 

Sinc e th e direct d etec tion of the magnetic field on 
8 1 OriC (Donati et al. 2002), significant magnetic fields have 


located near the Alfven radius ft : \, see equations 40 and 41 of 


ud-Doula et al. 


(2014). 


been measured on ~10 other O stars (see Petit et al. ||2013 


Morel et al. 2014 1 . Significantly, most of these stars are 
classified as Of?p (Walborn et al. 2010| ), though the rea¬ 
son for this connection between spectral morphology and 
magnetism is not yet fully understood. All magnetic O stars 
share physical and magnetic properties with 8 1 Ori C. As dis¬ 
cussed by Petit et al. (20131 all but one of the magnetic O- 
type stars are slow rotators, consistent with strong magnetic 


braking associated with their powerful stellar winds (ud- 


Doula et al. 20091. Their observed optical magnetospheric 


properties are well explained within the dynamical magne¬ 
tosphere framework: with negligible rotational support, the 
plasma environment in the magnetosphere should be quite 
dynamic, with trapped plasma falling back down onto the 
star in not much longer than a free-fall time, as confirmed 
by numerical MHD simulations (ud-Doula et al. 2013). 

Magnetic O stars are generally an order of magnitude 
more luminous in X-rays than non-magnetic O stars of 
the same spectral subtype, and their X-rays are somewhat 
harder, although the X-ray emission from other Of?p stars 


is generally not as hard as that of 6 )1 OriC (Naze et al. 
[2014). Given the similarity in physical and magnetic prop¬ 
erties among the magnetic O-type stars, this difference in 
spectral X-ray properties has yet to be understood. 

The star that is the subject of this paper, the slowly 
rotating O-star NGC 1624-2, has a field strength more than 
an order of magnitude larger than any other known mag¬ 
netic O-stars, at nearly 20kG (Wade et al.|2012a hereafter 


W12al. Given that field strength, and a similar expected 


wind mass-loss rate as 8 1 OriC and the Of?p stars, it should 
have a significantly larger magnetic closure radius and thus 
a significantly larger dynamical magnetosphere, extending 
out to ~ 11 ¥U. In order to explore the X-ray generation 
and also the possible magnetospheric absorption of X-rays in 
this star’s giant magnetosphere, we have obtained two Chan¬ 
dra ACIS observations of NGC 1624-2 at rotational phases 
that correspond to two viewing angles, respectively when the 
magnetic pole and equator are closest to the line of sight. 

In this paper we present detailed analysis of the X-ray 
spectral properties at both phases, including modelling of 
any possible local absorption by partially ionised plasma 
in the magnetosphere. This allows us to draw conclusions 
about the physical properties of the magnetically confined 
wind shocks, the wind feeding rate and heating efficiency, 
the post-shock cooling plasma, and the spatial distribution 
of magnetospheric material. In addition to this empirical 
spectral modelling, we present a comparison with recent 


MHD models of X-ray emission from magnetospheres (ud- 
Doula et al.||~2014l. Also preliminary results from the post¬ 


processing of MHD simulations of a generic O star magne¬ 
tosphere provide information about the expected viewing- 
angle dependence of the X-ray absorption column density 
(Owocki et al. in preparation). The goal of this paper is to 
constrain models and extend our understanding of the dy¬ 
namical magnetospheres of O stars via observations of the 
most extreme example currently known. 


2 STELLAR PROPERTIES OF NGC 1624-2 

NGC 1624-2 is an 07 f?cp star that was recently singled out 
as the most extreme example of the Of?p spectroscopic cat- 


© 2014 RAS, MNRAS 000, [TlfTs] 
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egory (Walborn et al. 20101. Soon after, a ~20kG field was 
derived by the MiMeS Collaboration (W12a). Although very 


strong magnetic fields have been measured for a handful 
of lower mass ApBp stars (e.g. 34 kG Babcock’s star; see 
Mathys et al. 1997), NGC 1624-2 has by far the strongest 


magnetic held ever measured on an O-type star: all the oth¬ 
ers have -Bpde < 2.5 kG. According to the properties de¬ 
rived by |W12a| summarised here in Tablc[l] NGC 1624-2’s 
magnetosphere is estimated to extend to Ra ~ lli?* hence 
trapping 95% of the outflowing wind, much more than other 
magnetic O stars which have Ra of just a few stellar radii. 


|W12a| determined a rotational period of 157.99 d from 
the observed spectral variations. This slow rotation is con¬ 
sistent with the circular polarisation signal variation, and 
with the narrow photospheric line profiles. As pointed out 


by Sundqvist et al. (20131, the width of the spectral lines, 


which exhibit Zeeman splitting, is completely accounted for 
by Zeeman broadening, and does not show any signs of 
the large macroturbulent broadening generally observed in 
both normal and magnetic O-type stars. This suggests that 
NGC 1624-2 is the only known magnetic O-type star with a 
held strong enough to stabilise its deep photosphere, provid¬ 
ing the first empirical constraints on the connection between 
the sub-surface convection zone present in models of hot, 


massive stars (Cantiello et al. 20091, and macro-turbulent 
line broadening. 

The dense magnetospheres of magnetic O stars pro¬ 
duce emission in optical recombination lines, such as the 
Balmer lines, which is modulated according to the rotational 
period. This variation is due to changes in the projected 
distribution of optically thick magnetospheric material (see 
Fig.|T| [Sundqvist et al. 20121. Although these slowly rotating 
magnetospheres are quite dynamic in nature, their volume- 
integrated properties are quite stable (Howarth et al.|[2007 


ud-Doula et al.|2013l from one rotation cycle to another. 


Although the detailed magnetic geometry of NGC 1624- 
2 has yet to be established, the single-wave, sinusoidal nature 
of the spectral line variation implies that we only view one 
magnetic hemisphere of a mainly dipolar field during the 


rotation of the star (W12a). However, the large variation in 
the Ha emission (Fig.|2j) between a “high state”, when the 
magnetic pole is most closely aligned with the observer’s 
line of sight, and a “low state”, when the magnetic equator 


is closest to the line of sight (e.g. Sundqvist et al. 20121 
indicates a large change in our viewing angle (see FigTjT]). 

NGC 1624 was observed for 10 ks with the Chandra X- 
ray Observatory in 2006 (ObsID 7473, PI Gordon). Even 
though NGC 1624-2 was detected with only 38 counts, only 
one count has an energy below 1 keV indicating rather hard 
emission. According to the analysis of |W12a| this spectrum 
could be reproduced with a hot (~ 2.3 keV) plasma, similar 
to that of 6 1 Ori C. It could also be reproduced with a cooler 
(~ 0.7 keV) plasma, similar to that of the other Of?p stars, 
that would be more heavily absorbed (Ah ~ 2 x 10 22 cm -2 ). 
We note that even during the low state, the Ha emission of 
NGC 1624-2 is much stronger than that of any other known 


magnetic O-type star (W12a see their Fig. 11). This sug¬ 
gests that the magnetosphere must indeed be very large, and 
higher S/N X-ray observations are the perfect tool to verify 
this hypothesis through the presence of X-ray absorption. 


Table 1 . Summary of adopted stellar, wind, magnetic, and mag¬ 
netospheric properties of NGC 1624-2 from |W12a| The stellar 
and magnetic parameters were derived from optical spectropo- 
larimetry. The unperturbed wind mass-loss rate M and terminal 
velocity Voo were calculated based on the theoretical prescriptions 
of ([Vink et al. 200l|. The wind magnetic confinement Alfven ra¬ 
dius i?Alf 5 an< l the Kepler co-rotation radius i?Kep are calculated 
as described by Petit et al.| ||2013||. 


Spectral type 

06.5f?cp - 08f?cp 

Teff (K) 

35 000 ± 2000 

log 9 (cgs) 

4.0 ± 0.2 

R* (R 0 ) 

10 ±3 

log (L*/L 0 ) 

5.10± 0.2 

M^ vo1 (Mq) 

28t£ 

log M b= o (M 0 yr -1 ) 

-6.8 

Voo,B=0 (kills -1 ) 

2875 

E(B - V) (mag) 

0.802 ±0.02 

R\ 

3.74 ± 0.1 

log d (pc) 

3.71 ± 0.1 

Prot (d) 

157.99 ± 0.94 

Pd (kG) 

~ 20 

Rar (P*) 

11 

R-Kep (P*) 

> 40 


3 OBSERVATIONS 


We obtained two observations of NGC 1624-2 with the bare 
ACIS-S camera in VFaint mode aboard the Chandra X- 
ray Observatory (PI: Petit). The ACIS-S camera allows for 
imaging spectroscopy combining moderate spectral resolu¬ 
tion (FWHMc^ 150 eV at 1.5 keV) and sub-arcsec spatial 
resolution. According to the ephemeris of W12a| (JD = 
[245 5967.0 ± 10] ± [157.99 ±0.94] x E), the observations were 
obtained at rotational phase 0.43 ±0.09 (ObsID 14572) and 
0.96 ±0.09 (ObsID 14571). These correspond roughly to the 
low state and the high state, respectively (see Fig[2|. 

The data were reprocessed according to standard reduc¬ 
tion procedure^ with ciao version 4.5. Using the cellde- 
tect tool, an X-ray source is detected at the position of 
NGC 1624-2 in both observations (Fig.[3|. The pile-up frac¬ 
tion, evaluated from maps of count rate per frame, is less 
than one percent for both observations. Table [2] compiles 
the aperture photometry calculated with the APRATES tool. 
The net photon flux is 40 percent higher during the Ha 
emission high state than at the low state. 


4 TIMING ANALYSIS 

For both observations, we produced light curves to search 
for variability on a time scale comparable or less than the 
50 ks exposure. Such a timescale corresponds to only 0.4 
percent of a rotation cycle, but is comparable to the wind 
flow timescale. The broad-band background-subtracted light 
curves, binned to an arbitrary 2 ks, are shown in red circles 
in Fig. [4] We applied the Gregory-Loredo variability algo¬ 
rithm to the source event list, as implemented in the Cl AO- 
4.5 routine GLVARY. No evidence of short term variability 
is found and the resulting probability weighted light curves 


2 The CIAO analysis threads are available at 
http://cxc.harvard.edu/ciao/threads/index.html 
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Figure 3. Chandra ACIS-S images of NGC 1624-2 (position marked by the green cross) during the Ha low state (left) and high state 
(right). The sources detected by CELLDETECT, including NGC 1624-2, are indicated by solid red ellipses. The spectral extraction source 
(thick) and background (thin) regions are shown as dashed blue circles. The other point sources circled in dashed blue were excluded 
from the background estimate. 



-0.5 0.0 0.5 1.0 1.5 

Rotational phase 


Figure 2. Equivalent width of the Ha line of NGC 1624-2 phased 
with the ephemeris of |W12a| The symbol shapes represent the 
rotation number since the ephemeris JDq. The optical spectra 
were obtained with ESPaDOnS at CFHT over a period of time 
corresponding to nearly 3 rotations of the star, by the MiMeS 
Collaboration ( [Wade et al. 20lT| and on PI time (Wade). The 
Chandra observations, indicated by the grey shaded area, were 
obtained during rotation cycle no. 3. 

are shown as black squares. We also applied a maximum- 
likelihood block algorithm that divides the event list into 
blocks of constant count rate (Wolk et al. 2005|. Each epoch 
is reproduced by a single block (grey bar), which also in- 


Table 2. Observation details and aperture photometry results. 



Hq low state 

Ha high state 

ObsID 

14572 

14571 

Date 

2013-08-03 

2013-10-27 

Start HJD 

2456508.20 

2456592.84 

Rotational phase 

0.43T0.09 

0.96T0.09 

Net exposure (ks) 

49.5 

49.5 

Net counts 

301.^29 

4 OI +35 

4cU -35 

Count rate (cts —1 ) 

6 .l 3 ;| x 10“ 3 

00 

1 + 
0 0 

X 

O 

1 

w 

Photon flux (Phcm - 2 s -1 ) 

1- 7to".i x 10-5 

2.43:2 x ! 0 -5 


dicates a lack of short-term variability. A comparison of 
the two epochs however confirms the statistically significant 
variability on the longer rotational timescale. 


5 SPECTRAL ANALYSIS 

We extracted the spectra of NGC 1624-2 and associated in¬ 
strumental responses using the standard SPECEXTRACT pro¬ 
cedure. We grouped the spectral channels to obtain a signal- 
to-noise ratio of 3 in each bin, allowing the use of % 2 statistics 
in our model fitting. Fig. [ 5 ] shows the resulting spectra (top 
panel) for the low state (thin black) and high state (thick 
red) observations. The bottom panel shows a representa¬ 
tion of the spectral response functions for both observations. 
They are practically identical, allowing meaningful direct vi¬ 
sual comparison of the spectra. In both observations, most 
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Low state 



High state 



Figure 4. Light curves between 0.3 and 10 keV, binned by 2ks (red circles). The probability-weighted light curves from the Gregory- 
Loredo variability algorithm implemented in the GLVARY routine are shown with black squares. The maximum-likelihood block algorithm 
resulted in a single block for each observation, shown by the grey bars. 



1 10 
Energy (keV) 


Figure 5. ACIS-S spectra of NGC 1624-2 during the low state 
(thin black) and the high state (thick red). The bottom panel 
shows a representation of the instrumental response, i.e. the spec¬ 
tra that would be observed if the emission model was flat. The 
small differences between the two epochs are caused by slight 
variations in response and adaptive signal-to-noise binning. 


of the photons have energies higher than 0.8 keV. The high 
state observation shows an excess of soft (<~1 keV) photons. 
Prominent spectral features, similar in both epochs, are seen 
between 1 and 2keV. 

We modeled the spectra with xspec 12.7 ( |Arnau d1996) 
with an optically thin plasma model from the Astrophysical 


Plasma Emission Code (APEC with atomdb 2.0.1, Smith 
|et al.|[2001| |Foster et al.|[20l2| ), which assumes collisional- 
radiative equilibrium (‘coronal emission’). Given the lack of 


detailed abundance analysis in the optical (W12al, we here 
adopted solar abundance values ( Asplund et al.|2009 l. 

X-rays in O-type stars are generally produced by 
shock-heated wind material that mostly cools by radia¬ 
tive emission. If the wind material is continuously heated 
and cooled, a time-averaged description of the wind should 
contain a continuous distribution of plasma temperatures 
(e.g. Antokhin et al. 2004). We therefore mimic such a 
distribution by using 4 plasma temperature components 
[T T T T 2) with fixed kT of 0.2, 0.6, 1.0, and 4.0 
keV, corresponding loosely to the peak emissivity of dom¬ 
inant spectral features likely to be present in the ACIS-S 
bandpass (H-like and He-like atoms of Ne, Mg, Si, S, and 
Fe). Given that these emissivity functions generally extend 
over a broad range of temperatures, we believe this approach 
is sufficient to model the characteristics of the emission spec¬ 
trum while minimising the number of fitting parameters (see 
discussion in Naze et al.|20 141. It is important to note that 
other representations of the post-shock plasma (e.g. different 
combinations of temperatures and their associated emission 
measures) can result in a similar emission spectrum, espe¬ 
cially at the relatively low resolution of ACIS-S. However 
the goal of this analysis is to investigate the global proper¬ 
ties of the phase-dependent absorption, and not the detailed 
distribution of emission measure. 


We attenuate the total emission model spectrum de¬ 
scribed above with the cold ISM absorption model tbabs 
[ColdAbsx (4T)] of Wilms et al. (20001. We adopt an ISM 
column density IVjj = 0.48 x 10^ cm -2 , determined from 
E(B — V) = 0.804, Ry = 3.7 4 (|W12a|), which results in 
Ay = 3.01 |Giiver fe 6zel||2009| ). 


In addition to this “ISM only” model, we consider the 
possibility of extra local circumstellar absorption (“LOC”) 
in two ways. First, we consider an additional cold absorption 
(“Cold LOC”) with the same absorption model as the ISM 
[ColdAbsx (ColdAbsx (4T))]. Second, we consider the pos¬ 
sibility that this additional absorber is composed of warm 
wind material at near the ~35kK effective temperature of 
the star (“Warm LOC”; [ColdAbsx ( WarmAbsx (47 1 ))]). Be- 
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cause H and He are mostly ionized at this temperature, the 
wavelength-dependent opacity will be different than for a 
cold absorber model, transmitting more flux at lower ener¬ 
gies. We therefore use the model slabtabs provided in the 
contributed xspec package windprofili]^] Slabtabs models 
absorption by a simple slab of ionized material, with elemen¬ 
tal abundance and ionization balance chosen for the appli¬ 


cation, and using opacities from Verner & Yakovlev (19951 


as implemented by Leutenegger et al. (20101. We use their 
opacity model calculated for a generic O-type star wind with 
solar abundances (see their Fig. 6) 

We consider a simplified geometric model where the 
absorber is located between the emission site and the ob¬ 
server. In this case, the attenuation goes as e ~ TX , where t\ 
is the wavelength-dependent optical depth. It is likely that 
the sites of X-ray emission are embedded within the warm 
wind material. Such geometric models exist for spherically 


symmetric non-magnetic winds (e.g. windtabs Leutenegger 
et al. 20101, however they might not be readily applicable 


here given that a magnetosphere is not expected to be spher¬ 
ically symmetric. 

We proceed to model the observations using two differ¬ 
ent approaches, which correspond to the results presented 
in the following two subsections: 

(i) In the simplest case (“individual models”), we let all 
model parameters (column density nH hOC and the distance- 
weighted emission measure EMi associated with each tem¬ 
perature component) vary independently from one observa¬ 
tion to the other. More specifically, this assumes that the 
plasma temperature distribution (the 4T components) can 
differ from one epoch to another. This could be due to intrin¬ 
sic changes in the temperature distribution with time, or by 
occultation of spatially stratified temperature distribution 
(see below). 

(ii) The spectral variations of magnetic stars are usually 
understood to be due to a change in the observer’s view¬ 
ing angle as the tilted magnetosphere co-rotates with the 
stellar surface. Although these magnetospheres are dynamic 
in nature, there are many indications that their volume- 
integrated properties are quite stable, namely based on nu¬ 
merical simulations ( ud-Doula et al.|2013 ), on observational 
diagnostics ( Howarth et al.|2007 Townsend et al. 20131, and 
on the lack of short term variability observed here. With 
this picture in mind, we therefore do not expect the intrin¬ 
sic X-ray emissivity of the magnetosphere itself to change 
significantly in time, but rather its sky-projected structure 
to modulate with the rotational phase. We thus also at¬ 
tempt to simultaneously model both epochs with a common 
4T emission model (“Joint models”). We allow for a change 
in the overall flux of the emission model spectrum during 
the low state [Constantx (AT)] to account for a possible oc¬ 
cultation of the circumstellar X-ray emitting plasma by the 
stellar disk when the magnetic field is seen near equator-on 


I Gagne et al. 2005). Since the column density depends on 
the ray’s path, and therefore on the viewing angle, the local 
absorption is also allowed to differ between the two epochs. 
This, however, makes the simplifying assumption that there 
is no strong spatial gradient of plasma temperature inside 


3 http://heasarc.nasa.gov/xanadu/xspec/models/windprof.html 


the magnetosphere, which would require more precise mod¬ 
els to constrain. In summary, the seven fitted parameters 
for the joint model are the column densities for the low and 
high states, the distance-weighted emission measure associ¬ 
ated with each temperature component, and a flux reduction 
factor / for the low state. 


5.1 Results from individual models 

We fit the individual models, with the three different absorp¬ 
tion models (“ISM only”, “Cold LOC”, and “Warm LOC”), 
to each observation independently. The top panels of Fig. [6] 
show the best fit model spectra folded through the ACIS- 
S response function (red, green and blue lines) overlaid on 
the observed spectra (grey histogram). The 3 bottom pan¬ 
els show, from top to bottom: (i) The attenuated model 
spectra as would be seen from Earth, with dominant spec¬ 
tral features identified, (ii) The model spectra corrected for 
the ISM absorption, as emergent from the magnetosphere, 

(iii) The emission model spectra (AT components; essen¬ 
tially both ISM and LOC corrected) associated with these 
best fit models. 

The fit results are compiled in Table[3] Columns 1-3 de¬ 
scribe the absorption model under consideration. Columns 
4-7 give the normalisation of each T component, and column 
8 gives the mean plasma temperature weighted by this nor¬ 
malisation. The energy fluxes in column 9 are calculated by 
integrating the model spectra between 0.3 and 10.0 keV. We 
also parametrize the relative hardness of the model spectra 
by measuring a flux hardness ratio (column 10) that com¬ 
pares soft (0.3-2.0 keV) and hard (2.0-10.0 keV) bands. The 
reduced x 2 of the model fits are given in column 11. We also 
provide the reduced \ 2 values that correspond to the la, 
2a, and 3a confidence intervals for the degrees of freedom 
of each observation (23 and 33 for low state and high state, 
respectively). 

In the second part of Table [3j we compile the character¬ 
istics of the the ISM-corrected model spectra (columns 12- 
15), and those of the emission 4T model spectra (columns 
16-19). In addition to the flux and hardness ratio described 
above, we also calculate the associated luminosity in the X- 
ray band log(Lx) and the ratio of the X-ray luminosity to 
the bolometric luminosity log(Z/x/Tboi). 

From these models, we reach the following conclusions: 

(i) For both observations, the models without local absorp¬ 
tion (ISM only) provide a poorer fit than those with extra 
local absorption. 

In the case of the low state, the improvement in the fit 
is statistically significant. In the case of the high state, all 
three models are within the la confidence interval. Indeed, 
the emission spectra (bottom panels of Fig. [6| and fluxes 
(column 16 of Table [3| of the different models are more 
similar during the high state observation than during the 
low state observation. 

The need for this local absorption component can be un¬ 
derstood by the fact that a hot plasma component is neces¬ 
sary to explain the hard (>3keV) tail of the observed spec¬ 
tra. However, the model only absorbed by the ISM material 
is too hot (feT avg ~3.7 keV) to accurately reproduce the very 
prominent Mg and Si features, which have peak emissivi- 
ties at cooler temperatures (logT~6.8 and 7.0 or fcT~0.5 
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Figure 6. Individual fits to the low state (left) and the high state (right). The emission component of the model is a 4-component 
fixed-temperature apec spectrum (0.2, 0.4, 1.0 and 4.0 keV). The absorption models are a cold ISM only (thin red), a cold ISM plus 
cold local absorption (green), and a cold ISM plus warm local absorption (thick blue). The three lower panels show, from top to bottom, 
the best fit for each attenuation model, the best fit local absorption component without the ISM absorption, and the best fit emission 
component (apec only; no ISM/LOC absorption). Note the prominent Mg and Si spectral features between 1 and 2 keV, that can only 
be reproduced with a relatively cool plasma component. 


and 0.8keV, respectively). The inclusion of a cooler compo¬ 
nent with a stronger emission measure to model the Mg and 
Si features cannot reproduce the lack of flux at low energy 
(< 0.8 keV) without the presence of additional absorption. 

(ii) During both epochs, the amount of local absorption is 
very large. 

One explanation for this result could be that the ISM 
column density has been underestimated. However, the Ay 
corresponding to these additional column densities would 
be ll±3mag and 6±3mag for the low and high state, re¬ 
spectively. Even during the high state, this difference is 
larger than the uncertainties on the reddening determined 
bv IW12al (see Tablc[l]|. It is therefore rather unlikely that 
the extra absorption is of interstellar origin. 

The two models with cold and warm circumstellar ab¬ 
sorption result in statistically similar fits. If we consider the 
likely hypothesis that this absorption is the result of warm, 
ionised material trapped in the magnetosphere, the column 
densities needed to reproduce the observed spectra are quite 


large. As a base for comparison, although the expected wind 
geometry is different, such a column density would rival the 
characteristic column density for the non-magnetic early-0 


supergiants if Pup or HD 150136 (sec Table 1 of Lcutcncggcr 
et al.|2010 l. However, as pointed out by W12a NGC 1624-2 


has by far the largest Ha equivalent width of all the magnetic 
O-type stars, implying a greater amount of magnetospheric 
plasma. 


(iii) There are large differences in the emission model 
spectra when comparing the low and high state observations 
considering the same absorption model. 

The model spectra in the low state are softer than in the 
high state, as shown in Table [3] by the hardness ratio in col¬ 
umn 17 and the fcTkvg in column 8. The flux of the emission 
component is higher during the low state by a factor 3-5. The 
conclusion is the same when considering a cold or a warm 
circumstellar absorption. The emission spectra correspond¬ 
ing to a different absorption model for a given observation 
are only slightly different, as can be seen by comparing the 
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Table 3. Individual models (shown in Fig. |§|. The distance-weighted emission measure normalization of the APEC components ( EMi ) 
are in units of 10 -14 cm —5 . These should be multiplied by And 2 to be converted into volume emission measures in cm -3 . We here compile 
these factors with a greater accuracy than their typical formal modeling uncertainty, in order for our model spectra to be reproducible. 


The formal error on the Earth-observed flux is or the order of 1 x 

10 14 ergs 

2 . All other quantities are given to the significant 

digit. 











Model 

«H ism 

nff LOC 

EMi 

EMi 

EM$ 

EM4. 

fcT avg 

Flux 

H/S 

X r 2 ed 


10 22 cm- 2 

10 22 cm —2 

(0.2 keV) 

(0.6 keV) 

(1.0 keV) 

(4.0 keV) 

(keV) 

erg s -1 cm —2 



(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 



Low state % 2 ed (l cr ?2cr,3f7)= 

1.12, 1.55, 2.00 for 23 DF 




ISM only 

0.48 

0 

3.59E-06 

0 

0 

4.06E-05 

3.7 

4.7E-14 

2.5 

2.12 

Cold LOG 

0.48 

1.78T0.50 

6.83E-04 

4.28E-05 

6.14E-05 

3.18E-05 

0.4 

4.6E-14 

2.3 

1.52 

Warm LOC 

0.48 

3.93T0.73 

1.91E-03 

9.85E-09 

1.12E-04 

2.71E-05 

0.3 

4.6E-14 

2.0 

1.14 



High state x? ed (1 ct,2ct,3<t)= 

=1.10, 1.45, 

1.81 for 33 DF 




ISM only 

0.48 

0 

0 

0 

4.14E-06 

5.82E-05 

3.8 

7.0E-14 

2.3 

1.09 

Cold LOC 

0.48 

0.92T0.44 

1.78E-04 

2.33E-05 

3.65E-05 

5.61E-05 

1.1 

7.1E-14 

2.3 

0.82 

Warm LOC 

0.48 

1.96T0.46 

5.71E-05 

3.34E-05 

4.98E-05 

5.46E-05 

1.5 

7.1E-14 

2.2 

0.79 



ISM corrected 



Emission component only 



Model 

Flux 

H/S 

log (Lx) log(Lx/L bo i) 

Flux 

H/S 

log(L x ) 

log(Lx/L bo l) 




ergs -1 cm -2 


ergs 1 


ergs -1 cm' 

-2 

ergs -1 




(1) 

(12) 

(13) 

(14) 

(15) 

(16) 

(17) 

(18) 

(19) 







Low state 







ISM only 

7.0E-14 

1.0 

32.3 

-6.3 

7.0E-14 

1.0 

32.3 

-6.3 



Cold LOC 

5.5E-14 

1.5 

32.2 

-6.4 

1.1E-12 

0.04 

33.6 

-5.1 



Warm LOC 

7.3E-14 

0.8 

32.4 

-6.3 

2.6E-12 

0.02 

33.9 

-4.8 







High state 







ISM only 

1.0E-13 

1.0 

32.5 

-6.2 

1.0E-13 

1.0 

32.5 

-6.2 



Cold LOC 

9.0E-14 

1.3 

32.5 

-6.2 

4.4E-13 

0.1 

33.1 

-5.5 



Warm LOC 

9.6E-14 

1.1 

32.5 

-6.2 

3.4E-13 

0.2 

33.0 

-5.7 




green and blue curves in the bottom panels of Fig. [6] (also 
column 16 and 17 of Table[3|, the fluxes and the spectral 
hardness are different by a factor of < 2. 

Such a result would be difficult to explain in terms of oc- 
cultation of a spatially stratified temperature distribution. 
As discussed in the following section, simultaneously mod¬ 
elling both epochs with a common emission model provide a 
slightly poorer but still adequate fit to the observations. A 
high-resolution X-ray spectrum of NGC 1624-2 would pro¬ 
vide a useful constraint for the detailed properties of the 
magnetosphere emission processes. 

(iv) The X-ray emission of NGC 162f-2 is very luminous, 
and heavily attenuated. 

Depending on the characteristics of the additional ab¬ 
sorption component, the emergent (ISM corrected) X-ray 
luminosities between 0.3 and 10 keV are log(Lx) ~32.3 
and 32.5 for the low state and high state, corresponding to 
log(Lx/F B oi) of -6.3 and -6.2, respectively, similar to what 
is seen in other magnetic O-type stars (Naze et al. |2014[ ). 
These models also imply that the shock-heated material pro¬ 
duces high intrinsic X-ray fluxes, corresponding to roughly 
log(Lx) of ~ 33.7 and ~ 33.0, or log(Lx/I/Boi) of —5.0 and 
—5.6. This would imply that 70-95 percent of the MCWS- 
produced X-ray emission is absorbed before escaping the 
magnetosphere. As a comparison, |Cohen et al. (2011) de¬ 
termined that 80 percent of X-ray emission produced in the 
wind of the non-magnetic star HD 93129A, one of the earli¬ 
est known O stars with a mass-loss of nearly 10 _5 Mq yr _1 , 
is absorbed before escaping the wind. We consider the val¬ 


ues derived here to be first order estimates, as the geometry 
of the models might be too simplified to accurately describe 
the structure of a real magnetosphere. 


5.2 Results from joint modelling 

Let us now turn our attention to the joint models described 
previously, where the 4 T emissivity model is constrained to 
simultaneously fit both low and high state spectra, allow¬ 
ing only for a change in overall flux reduction during the 
low state to account for the occultation of the shock-heated 
plasma by the stellar disk. The results of these joint fits, 
for the three absorption models, are compiled in Table [4] 
The blank entries represent values that were kept common 
between the two epochs. 

Once again the models with extra local absorption bet¬ 
ter reproduce the observed spectra, in a statistically signifi¬ 
cant way for the low state (see column 11). The joint models 
provide marginally poorer fits than the individual models, 
although all the models including local absorption are within 
the 2a confidence interval. Fig.[7]shows a comparison of the 
best individual models and the joint model for the warm 
local absorption model. During the low state (left panel), 
the joint emission model spectrum (bottom, in thick blue) 
is harder than the best individual emission model (in thin 
red). During the high state (right panel), the joint emission 
model spectrum is softer than the best individual emission 
model (in thin green). This can also be seen by comparing 
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Table 4. Same as Table [ 3 ] for joint models (Fig. [ 7 ] shows a comparison between the joint and individual models for the warm local 


absorption case). 

Model 

u h ism 

nH LOC 

EM 1 

EM 2 

EMg 

EM 4 

fcT avg 

Flux 

H/S 

CL 


10 22 cm -2 

10 22 cm —2 

(0.2 keV) (0.6 keV) 

(1.0 keV) 

(4.0 keV) 

(keV) 

ergs -1 cm —2 



(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 



Low state 

X?ed (1 o-,2ct,3o> 

=1.12, 1.55, 2.00 for 23 DF 




ISM only 


0 






4.6E-14 

2.5 

2.28 

Cold LOC 


1.41T0.36 






4.9E-14 

2.6 

1.68 

Warm LOC 


3.02T0.40 






4.9E-14 

2.4 

1.35 



High state 

X? ed (l<r,2<7,3<x): 

= 1.10, 1.45, 

1.81 for 33 DF 




ISM only 

0.48 

0 

0 

0 

1.57E-06 

6.25E-05 

3.9 

7.3E-14 

2.496 

1.13 

Cold LOC 

0.48 

1.13T0.36 

3.57E-04 3.83E-05 

4.63E-05 

5.06E-05 

0.7 

6.8E-14 

2.077 

0.85 

Warm LOC 

0.48 

2.61T0.40 

4.39E-04 3.61E-05 

8.29E-05 

4.61E-05 

0.6 

6.8E-14 

1.994 

0.86 



ISM corrected 



Emission component only 



Model 

Flux 

H/S 

log (Tx) 

log(L x /hboi) 

Flux 

H/S 

log(Lx) 

log(Lx/T bo i) 




ergs -1 cm - 

-2 

ergs -1 


ergs -1 cm - 

-2 

ergs -1 




(1) 

(12) 

(13) 

(14) 

(15) 

(16) 

(17) 

(18) 

(19) 







Low state 







ISM only 

6.5E-14 

1.1 

32.3 

- 6.4 

0.6 Fh igh 


32.3 

-6.4 



Cold LOC 

5.9E-14 

1.7 

32.3 

- 6.4 

0.8 F high 


33.2 

-5.4 



Warm LOC 

6.7E-14 

1.1 

32.3 

- 6.4 

0.8 -Fh igh 


33.3 

-5.3 







High state 







ISM only 

1.0E-13 

1.1 

32.5 

- 6.2 

1.0E-13 

1.1 

32.5 

-6.2 



Cold LOC 

8.6E-14 

1.3 

32.4 

- 6.2 

7.1E-13 

0.08 

33.4 

-5.3 



Warm LOC 

1.0E-13 

0.9 

32.5 

- 6.2 

8.7E-13 

0.07 

33.4 

-5.2 




the hardness ratios of the emission models (column 17) and 
fcTkv g (column 8) of Table [ 3 ] and [ 4 ] 

Because of these differences in emission model spectra 
between the individual and joint approach, the joint model 
for the low state requires less absorption to reproduce the 
lack of soft flux (and slightly underestimates the flux in the 
Mg and Si features) compared to the individual models. The 
opposite effect is seen for the high state. This leads to a 
less pronounced difference in the column density of the local 
absorption between the two epochs (column 3). 

Given the similar quality of the fits from the individual 
and joint models, the emergent, ISM-corrected fluxes and 
luminosities are similar (columns 12 and 14). Considering 
now the models with local absorption, the intrinsic flux of 
the emission component during the high state (column 16) 
is higher than for the individual models by a factor ~2, 
leading to an intrinsic luminosity of log(Lx) ~ 33.4. With 
this simplified model, the best value for the flux reduction 
(Tniin ~ 0.8 F 1 n „) during the minimum phase implies that 
the stellar disk must therefore occult percent of the X- 
ray emitting material when the magnetosphere is seen close 
to edge-on. In addition, the circumstellar material attenu¬ 
ates 85-90% of the X-ray emission at both phases. 


6 DISCUSSION 

In the pioneering work of Babel fc Montmerle| (1 997bj ), 
the X-ray luminosity originating from magnetically confined 
wind shocks fundamentally scales linearly with wind mass- 
loss rate, because the thermal energy content of the radi¬ 
ating post-shock plasma depends on the wind kinetic en¬ 


ergy flux across the magnetically channeled shock front. This 
model was recently refined by ud-Doula et al. (2014), who 
combined MHD simulations and analytical scalings taking 
into account the effects of instrument-specific bandpass and 
“shock retreat”. 

Shock retreat occurs in low density environments, where 
the cooling length of the shocked material becomes compa¬ 
rable to the length of the magnetic loop. In this case, the 
X-ray emission is expected to be be spatially extended above 
and below the magnetic equator, which lowers the pre-shock 
velocity and thus the kinetic energy available to power the 
shocks. 

The overall importance of the shock-retreat effect can 
b e characterised by the c ooling parameter defined in eqn. 25 
of 


ud-Doula et al. 


(2014|f| 


: 0.034 


VgRl 2 
M-8 ’ 


(1) 


where vg = Uoo/(10 8 cms -1 ), R 12 = R*/{ 10 12 cm), and 
M- 6 = M/(10~ 6 M@ yr -1 ). The case where x > 1 corre¬ 
sponds to a regime where shock retreat is important. The 
case where x < 1> represents a regime with dense condi¬ 
tions where the cooling length is small compared to the size 
of the magnetic loop, and the X-ray emission remains from 
near the loop apex at the magnetic equator. 

Using the values listed in Table[I] we find x°o ~ 10 


4 It is important to note that the models of ud-Doula et al.]|2014) 
are parametrized by the mass-feeding rate at the base of an un¬ 
perturbed wind, which is not the same as the actual mass-loss 
ulimately lost by the star. 
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Figure 7. Joint fit of the minimum phase (left) and the maximum phase (right) using a common emission model for both epochs, 
composed of a 4-temperature apec model (with components at 0.2, 0.4, 1.0 and 4.0 keV). Only the overall normalisation of the emission 
spectrum is allowed to change between the epochs. The absorption model shown here is the ’’cold ISM with warm local absorption” 
model. The best-fit joint model is shown in thick blue, compared with the best fits to the observations taken individually from Fig. [6] 
(in thin red for the low state and thin green for the high state) for the same absorption model. The top panel shows the models folded 
through the instrumental response compared to the data, and the bottom panels shows the models, the models corrected for the ISM 
absorption, and the emission models (i.e. j T only). 


for NGC 1624-2, implying a moderately strong shock re¬ 
treat, with associated moderation of X-ray emission level 
and hardness. Given the strong dependence on the rather 
uncertain wind terminal speed, a slight reduction to a lower 
Vaa value could result in a cooling parameter as low as unity, 
but likely not much lower. Fig.[8] shows the expected spa¬ 
tial distribution of the integrated X-ray emission between 
0.3 and lOkeV from the analytical models described by |ud-| 


Doula et al. (2014), for the Xoo = 1 and — 10 cases (see 


also Owocki et al. in preparation). Although the latitudinal 
distribution of the X-ray emission is strongly dependent on 
Xoo, its radial extent is concentrated relatively close to the 
star (~ 1.5-3 R *), independent of Xoo as well as the closure 
radius, located near Ra . This can be understood as due to 
the competition between the increase of the velocity shock 
jump with distance from the stellar surface, and the density 
decrease due to the areal divergence of the field loop. 


Our joint modelling of the NGC 1624-2 spectra suggests 


that the stellar disk occults ~20 percent of the X-ray emit¬ 
ting material. Using a simple geometric model where the X- 
ray emission originates from a ring concentrated at a single 
radius Rx, the maximum occultation (when the magneto¬ 
sphere is perfectly edge-on) can be expressed as: 




7T 



( 2 ) 


Fig.0 shows the constraints on Rx from this simple occul¬ 
tation model, suggesting that the bulk of the X-ray emis¬ 
sion is indeed generated relatively close to the star (~2 R+), 
compared to the extent of the largest closed magnetic loops 
(~10 R+). 

In the idealized limit of 100% shock efficiency, the scal¬ 
ing analysis presented by |ud-Doula et al.| ( |2014| ) predicts a 
theoretical X-ray emission for NGC 1624-2 of log(Lx) = 34.4 
in the 0.3 - lOkeV bandpass. By comparison, the empirical 
analysis presented here gives a factor ~ 10 lower luminos¬ 
ity for the intrinsic emission component in this bandpass 
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Figure 8. Relative spatial distribution of the X-ray emission 
according to the analytical models of |ud-Doula et al.] ( |2014| see 
their appendix A) for cooling parameters Xoo = 1 (top) and Xoo = 
10 (bottom). The shock retreat distributes the X-ray emission 
over a wider range of latitude. However, the bulk of the X-ray 
emission is always located relatively close to the star, even for 
large closure radius (here 4it*). 


(e.g. log Lx = 33.4 for the joint model with warm local ab¬ 
sorption). This indicates a shock efficiency factor of ~10%, 
in good agreement with the estimates obtained from com¬ 
parison with MHD simulations (ud-Doul a et al.| 2014l, and 
also suggested by observations of magnetic massive stars 
(Naze et al. 20141. This supports our conclusion that the 
X-ray emission from NGC 1624-2 is highly attenuated, with 
an amount similar to what is seen in the winds of much 
earlier non-magnetic O-type stars. 

The high intrinsic X-ray luminosity of NGC 1624-2 is 
therefore consistent with the strong X-ray emission seen in 
other magnetic O-type stars. Empirically, the observed spec¬ 
trum of NGC 1624-2 is hard like that of 0 1 OriC; but intrin¬ 
sically it is softer, making it more similar to the other Of?p 
stars. On the other hand, NGC 1624-2 is distinguished from 
the other Of?p stars by its large (phase-dependent) absorp¬ 
tion, due to its spatially large magnetosphere. 

Direct MHD simulations for high degrees of magnetic 
wind confinement, as are needed for NGC 1624-2, are not 
computationally feasible at this time, although new de¬ 
velopments in rigid-field hydrodynamics models might be 


available in a near future (Townsend et al. 2007 Bard & 
|Towns cnd 2015]). However, we can use the 3-D MHD simula¬ 
tion presented by ud-Doula et al. (20131 for a lower magnetic 



R* / R. 


Figure 9. Reduction of the observed emission (/ = F 0 t, a /Ftot) 
caused by the occultation by the stellar disk of a ring of X-rays 
located at a radius Rx (black thick curve). The thin middle line 
and shaded area show the constraints obtained by the joint mod¬ 
elling with the “warm” circumstellar absorption. This supports 
the notion that the bulk of the X-rays are emitted relatively close 
to the stellar surface, as compared to the radius of the outermost 
closed loops (see text for details). 


wind confinement, similar to that found for other magnetic 
O-type stars, as a baseline to explore the level and modula¬ 
tion amplitude of X-ray absorbing column densities. 

Fig. |10| shows contour plots of the column density Nh, 
time-averaged over 50 ks, between any location in the mag¬ 
netic equatorial plane and the observer. The right panel 
corresponds to the high state, when the magnetosphere is 
viewed pole-on. The observer is therefore located in front of 
the page. The left panel corresponds to the low state, when 
the magnetosphere is viewed edge-on - the observer is there¬ 
fore located at the top of the page and viewing in the plane of 
the page toward the bottom. The grey region at the bottom 
marks the occultation of a portion of the equatorial plane 
by the stellar disk. The green dashed circles mark 2 L*, the 
expected location of the bulk of the X-ray emission. 

This 3D simulation predicts a column density between 
0.1-0.5X 10 22 cm' 2 between the observer and X-ray emission 
originating at ~ 2 i?*. This is consistent with the low ab¬ 
sorption column observed in other magnetic O-type stars 
(e.g. Gagne et al.|2005 Naze et al.|2014 ). Regions of higher 
column density exist along the green circle for the low state 
on the far side of the magnetic equatorial plane, but these re¬ 
gions are already occulted by the stellar disk. It is therefore 
likely that no large changes in absorption will be observed 
for stars with magnetospheres significantly smaller than that 
of NGC 1624-2. 

The empirical column density for NGC 1624-2 is higher 
(l-4x 10 22 cm -2 ), which, in these simulations, only occurs 
for rays that impact very close to the star. A suitably larger 
magnetosphere could produce such a larger attenuation, as 
a larger fraction of the wind will be channeled into closed 
magnetic loops, increasing the overall column density. An 
extension of the X-ray scaling analysis of |ud-Doula et ah] 
12014) to account for the spatial distribution of pre-shock 
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Figure 10. Time-averaged (50 ks) column density njj between any location on the magnetic equatorial plane and the observer, calculated 
from the 3-D MHD simulation presented by[ud -~Doula et al.| |2013]| for a low magnetic wind confinement, similar to that found for other 
magnetic O-type stars. In both panel, the magnetic equatorial plane is in the plane of the page, and the magnetic pole is out of the page. 
The right panel correspond to the high state when the magnetosphere is viewed pole-on (the observer is located in front of the page), 
and the left panel correspond to the low state when the magnetosphere is viewed edge-on (observer is located at the top of the page, 
looking into the page). The grey circle is the stellar disk and the grey zone at the bottom of the right panel marks regions occulted by 
the stellar disk. The dashed green circle is located at 2 R t . 


and cool post-shock material could help quantify the X-ray 
absorption for large magnetospheres, currently out of reach 
of present MHD simulations (Owocki et al. in prep). 


7 CONCLUSION 

We observed NGC 1624-2, the O-type star with the largest 
known magnetic field ( B p ~ 20 kG) using the ACIS-S cam¬ 
era on-board the Chandra X-ray Observatory. Our two ob¬ 
servations were obtained at the low and high state of the 
periodic Ha emission cycle, corresponding to the rotational 
phases where the magnetic field is the closest to equator-on 
and pole-on, respectively. 

In both observations most of the photon counts have 
energies higher than 0.8 keV. No evidence of short term 
variability is found within each exposure. However, the net 
photon flux is 40 percent higher during the Ha emission 
high state, which shows an excess of softer (~1 keV) photons 
compared to the low state. 

We model the spectra with optically thin, collisional ra¬ 
diative plasma models from the Astrophysical Plasma Emis¬ 
sion Code. We find that for both observations a large local 
magnetospheric absorption, in addition to that from the in¬ 
terstellar medium, provides a much better fit to the features 
present in the spectra. 

The X-ray emission of NGC 1624-2 is therefore very 
luminous, and heavily attenuated, by ~ 70 — 95 percent. 
Therefore, although the spectrum of NGC 1624-2 appears 
hard like that of the archetypical magnetic star 8 1 OriC, it 


is in fact intrinsically softer and more similar to the other 
magnetic Of?p stars. However, the large amount of circuin- 
stellar absorption present in NGC 1624-2 is not seen in any 
of the other magnetic stars. We argue that such a high at¬ 
tenuation, and its observed variation with stellar rotation, is 
consistent with a larger magnetic closure radius, within the 
framework of a dynamical magnetosphere in the Magneti¬ 
cally Confined Wind Shock X-ray context, and consistent 
with its remarkably strong magnetic field. 

The case of NGC 1624-2 also revealed that the plasma 
postshock temperature, and therefore the intrinsic X-ray 
spectral hardness, should not depend on the confinement 
level. The question thus remains as to why the emission of 
the archetype magnetic O-type star 9 1 Ori C is harder than 
that of the other magnetic O-type stars observed in X-rays. 
Upcoming HST/COS observations will address the ionisa¬ 
tion structure and kinematics of the plasma, as well as the 
feeding and mass-loss rates of the giant magnetosphere of 
NGC 1624-2. These new observations may reveal important 
information about the pre-shock wind velocity, that can be 
compared to 8 1 Ori C to verify whether the velocity struc¬ 
ture of the magnetosphere is responsible for this difference 
is plasma temperature. 
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